Quantum oscillations in the spin-boson model: 
Reduced visibility from non-Markovian effects and 
initial entanglement 

F.K. Wilhelm 

Institute for Quantum Computing and Department of Physics and Astronomy, 
University of Waterloo, 200 University Avenue W, Waterloo, ON, N2L 3G1, Canada 

E-mail: fwilhelm@iqc.ca 

Abstract. The loss of coherence of quantum oscillations is of fundamental interest 
as well as of practical importance in quantum computing. In solid-state experiments 
the oscillations show, next to the familiar exponential decay on time scales Ti/ 2 , an 
overall loss of amplitude. We solve the spin-Boson for a large class of initial conditions 
without the Markov approximation at the pure dephasing point. It is shown that a loss 
of visibility occurs in the form of a fast initial drop for factorized initial conditions and 
an overall reduction for entangled initial conditions. This loss of amplitude is distict 
from T2-decoherence with the difference being most drastic for environments with real 
or pseudo-gaps. This result is explained by bandwith effects in quantum noise as well 
as in terms of higher-order phase-breaking processes. For several experiments, such 
gapped environments arc identified. We conirm that this physics is valid beyond the 
pure dephasing point. 
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1. Introduction 



Quantum coherence driven by the goal of quantum computation is a central theme of 
present-day research in mesoscopic condensed-matter physics. In particular in the field 
of superconducting qubits [H El EH E] , spectacular successes have been achieved, such 
as Rabi oscillations, charge and flux echo, and a controlled-not gate [3J EJ U\ El M HH] • 
In order to achieve this, the decoherence due to the ubiquitous environmental degrees 
of freedom in solid-state systems had to be overcome. In fact, careful modelling of the 
environment allows to make predictions of the relaxation and dephasing times, T\ and 
T 2 , which are extracted either from spectroscopic line widths or from fitting exponential 
envelopes to Rabi- or Ramsey oscillation data [11] [12] . Theoretical predictions of T\ are 
in reasonable agreement with experiments. T 2 is mostly attributed to 1/ f- noise [SI IT2"] 
for which self-consistent theories can be formulated [T2l IT3] even though the detailed 
origins of that noise are not quantitatively predictable yet. Technically, the notion of 
exponential decay of coherences on scales Ti/2 is based on a Markov approximation of 
some, potentially implicit, type. 

On the other hand, a variety [HI El H3] of controlled experiments show an additional 
loss of visibility: Next to the exponential decay given by Ti/2, the amplitude of coherent 
oscillations is reduced even further, i.e., it does not extrapolate back to the full expected 
amplitude at t = 0. This behavior is not explained by a simple Ti/2-picture. This 
reduced visibility is an obstacle for the demonstration of macroscopic quantum effects 
in these devices, e.g., of the violation of Bell's inequality in a coupled qubit system 
[13] . and for quantum computing applications. It has so far mostly been attributed to 
technological shortcomings of the detector. This paper introduces a different, additional 
mechanism for a reduced visibility based on non-Markovian effects in the spin boson 
model. This generic effect is induced i) by higher order processes involving virtual 
intermediate states, which rapidly entangle system and environment and ii) potential 
initial entanglement between system and environment. We show, that the reduction 
originates from the off-resonant high-frequency parts of the environmental spectrum 
which do not contribute to Ti/ 2 . We demonstrate that reduced visibility is compatible 
with large T 2 -values in superohmic and gapped environments and has to be considered 
as an independent quantifier of decoherence. We identify such environments in recent 
experiments and estimate the loss of visibility they induce. 



2. The model 



As a generic model, we consider the spin-boson Hamiltonian [IB] 

H e = ~ (cos 9a z + sin 9a x ) + ^a z K(a* + 4) + ( a ! a * + ?>) • C 1 ) 

where the <jj are Pauli matrices and a and a) are Boson annihilation and creation 
operators. The bath is characterized by the spectral density J{uj) = J2i | | 2 <5(co? — Uj) 
which is related to the equilibrium spectral noise power S (uj) = J{oj) coth(u;/2T). This 
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model is realized in superconducting qubits, where the bath is the electromagnetic 
environment [TT] IT8] or by phonons [19] , which also play an important role in quantum 
dots (20]. Here and henceforth, we chose % = 1 and ks = 1- 

The spin-boson model is in general not exactly solvable. It has been treated with a 
number of approaches [TH [2TJ [22] . For quantum computing, the A, are small at u>i ~ E 
by design and the system-bath coupling is usually treated perturbatively. If the system- 
bath interaction also defines the longest time in the problem, a Markov approximation 
is justified (see, e.g., [18] for a recent review). This procedure leads to variants of the 
well-known Bloch-Redfield master equation which predicts strictly exponential decay of 
the spin projections contained in the system+bath density matrix p, Si = Tr [(<7j (g> l)p\ 
with time scales Ti/ 2 . Thus, conceivably, such an approach cannot describe the loss of 
visibility. Moreover, many designed environments use the option to allow for larger Aj 
at high frequencies [6], Ui^> E. 

In order to go beyond Bloch-Redfield, we use two approaches: For 9 = 0, the spin- 
boson model reduces to the exactly solvable independent boson model [231 EH ESI EE] . For 
9 7^ we will use perturbation theory in the qubit Hamiltonian to obtain approximate 
insights. We will recover similar physics in both cases. 



3. Pure dephasing point, 9 = 

We now proceed to the exact solution for a very general initial state in the independent 
boson limit. We can perform a Schmid decomposition of the initial density matrix in 
the qubit ® bath Hilbert space as 



„= ™ ™ . (2) 




As the bath is composed of noninteracting oscillators, the submatrices p^ = Y\ k p^' 
remain factorized within the bath modes k. Note that the submatrices do not have to 
be valid density matrices, only the complete p has to be. All matrices with bounded 
trace can be parameterized with the characteristic function \ °f the Wigner function 
X f {a<- k \aM*) = Tr[ ft ( ^W(-a< fe ))] with D^(a^) = ex P (a^a^ {k) - a* (fc) af } ) the 
displacement operator for mode k. Preparing such an initial state requires to use controls 
with 9{t) ^ at t < 0. 

The fact that for 9 = the Hamiltonians at different times t > commute 
allows us to exactly compute the propagator in interaction representation. The result 
can be written as a product of -D-operators, U(t) = Ilfc-D (/*fc(0) with Pk{t) = 
~zrSt ( elu}kt ~ 1)- The application of these to the density matrix is straightforward, the 
essentially induce a coordinate transformation on the a^: The matrix of characteristic 
functions will be 

= Tr fc (p^t = 0)D^)D^{-af)D^ k \pf)) (3) 
= e^Tr fc ( PlJ (t = 0)D^{-c^ + ^ - pf)) 
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where the phase factor comes from multiplying the displacement operators and 
reads 4$ = Im Piaff + afj 1 p* — p^p* . We can conclude that Xij\ a i a *it) = 

e^xf (o^\t),aff\t),0) without) = o^ + 5o^(t) with*4*>(f) = pf\t)-pf\t). 
This simple property that time evolution is a mere change of coordinate frame is the 
main rationale for resorting to x as a phase-space parameterization of the density matrix. 
The other rationale is that x makes it in particular easy to compute expectation values 
of qubit operators q using the corrolary Tr [pij] = Tr [p^- H k d\j\o) = l\ k X^(0, 0, t) as 

(g®i)=Tr[(g®i)p]= £ (i\q\j) H X ^(Sa(t), -Sa*(t), 0) (4) 

i,j£{0,l} k 

without any further integration. 

From knowing the full density matrix at any time t > we can calculate any 
property of the system we like, including correlation functions. The main purpose of 
this paper is, however, to discuss the quantum system alone in order to establish the 
connection to master equation approaches. 

Being at the pure dephasing point we can see that there is no dynamics on the 
diagonal elements: aa(t) = au(0) and (pa = 0. 

On the off-diagonal, things get more involved. We focus on studying the coherence 
given in terms of the charateristic Wigner function as 

s+ = ^{s x + is y ) = Y[xoi(Sa i(t),-5a* 01 (t),0). (5) 

k 

We can work out explicitly 

= (6) 

and it is easy to show that the phase factor drops out. 

We now apply this technique to a physical realistic realization by assuming a 
specific Xoi(<x, a*, 0). The main restriction on this function is that p has to be a valid 
density matrix, i.e. Hermitian, normalized, and positive. We are restricting ourselves 
to symmetrically entangling a qubit that is initially in an eigenstate of a x with classical 
states of the bath oscillators, i.e. by choosing for the complete initial density matrix 

$ = {uS%p$! (uK) (7) 

11 \ '33 

with p t h is the thermal density matrix for the bath and Uq — D^ k \S k ^a z ) displaces 
each bath mode in opposite directions in phase space conditionioned on the two states 
of the qubit, i.e. it performs a controlled unitary displacement by some mode-specific 
amound of the bath oscillators. Reusing the multiplication property of displacement 
operators this means that 

X g } («, a*, 0) = - 2S k \ a - 2d<*>) (8) 

with the phase space function of a thermal state being 

= e ~2aa* coth(u; fe /2T) ^ 
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K{t) = - - / -^S(cu) + l) 2 + v 2 {uj) + 1- (11) 



In this limit, we thus obtain 

s+ _ J-j- e -2coth(^ K /2T)|2d fe +M fe (t)| 2 _ ^q) 
fc 

taking the continuum limit and substituting c4 = (\k/2huk)(uk + ivk) with real 
dimensionless coefficients Ui and i>j. Taking the continuum limit we find = Res + = 
e~ K w cos e £ with 

1 Z" 00 dui 

2 Jo up' 
—2 ((1 + u{uj)) cos cut + smut)] . 

Different choices of u(uS) and v{uf) can lead to rich dynamics. We single out prominent 
quantities: The initial amplitude e~ K ^ = exp — / ^S(uj)(u 2 (uj) + v 2 (uj)) , the time 
constant Ti of exponential decay in the long-time limit, I/T2 = linif_ >00 d%K{€) = 
5(0) (u(0) + 1) and the visibility, the amplitude found when extrapolating the long-time 
behavior back to t — *■ 

= lim e^~ K{t) = e -^J^%)«WVH), (12) 

We analyze this result in important limiting cases. 
3.1. Factorizing initial conditions 

Factorizing initial conditions imply for the density matrix p(t = 0) = ps <8> Pb, i-e. in 
our notation tt(o;) = f (cu) = 0. This case corresponds to environments which are part of 
a detector or any piece of electromagnetic environment that is switched on during the 
experiment, e.g. the DC-SQUID detecting a flux qubit [TTJ, [271 12H] • It also applies to 
the case of excitons j29j [25] as the two-state system that are created by ultrafast laser 
control in the beginning of the experiment. From the above result, we can see that this 
is the only case in which the initial amplitude is unity, K(0) = 0, i.e. when the actual 
oscillations starts at full amplitude. In this case, K(t) can be written as |2H 130] 

roc 1 — ms(j;/ t^ r°° sin 2 \^r) 

K ®=L to-zpw-ih iM m - (13) 

The long-term expansion of this result can be written as K(t) = t/T2+logV[0,0]+O(l/t) 
with 1/T 2 = 5(0). Note, that the rate 1/T 2 is identical to the Bloch-Redfield result, but 
additionally the constant term V[0, 0] = exp — V / °° describes an overall loss of 
amplitude. Here V denotes the Cauchy mean value. 



The dynamics resulting from Eq. (13) can be interpreted in terms of environmental 



quantum noise. This is accomplished by looking at the second equality of eq. (13) 
which multiplies S[uS) with a spectral weight of width 1/t. When the system has been 
coupled to the environment for a time t, it samples its spectrum over a bandwidth of 
5uj = 1/t due to frequency-time uncertainty. This bandwidth goes to zero only in the 
long-time limit. Quantum-mechanically, the long-time limit only captures direct energy- 
conserving processes between system and environment, whereas at shorter times also 
higher-order processes involving energetically forbidden intermediate states play a role. 
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The leading process beyond T2 is the excitation of an environmental mode followed by 
the relaxation of another, excited state in the environment with infmitesimally different 
energy. These modes have to be distinct in order to leave a trace in the environment 
which is a necessary condition to entangle system and environment and hence cause true 
dephasing. Thus, it is crucial that the environment is in an excited state in the beginning 
of this cycle. In our case, this is guaranteed by the factorized initial conditions: The 
ground state of the spin-Boson Hamiltonian [16J is an appropriate dressed state of the 
spin. When the interaction is switched on, energy is redistributed: Forming this dressed 
state formes lowers the extra energy compared to the initial state. The extra energy is 
accomodated in bath excitations necessary for true dephasing as just explained. 

We will now study the consequences of these results for a number of important 
spectra, starting with the ones where T 2 — > 00, i.e. a constant amplitude at long times. 
These structured environments are an important test-bed for our approach because all 
the decoherence is from non-Markovian effects. The gapped Ohmic model approximately 

1 
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Figure 1. Dynamics of S x and its envelope in the gapped Ohmic model using a = 0.1 
lo c = 10E, ojir = 2E, and T = 0. The exact result is compared to the envelope which 
would be obtained within Bloch-Redfield theory for the same model. 

describes the effect of the quasiparticle transport channel shunting a Josephson junction 
[5T] . Its spectral density reads J g (u>) = aiUje^ U) ^ u)c Q{\uj\ — E g ) where is the Heaviside 
unit step function. In the limit of T <C E g <C u c , we find V[0, 0] = ( -^f-J e~" 17 where 
7 ~ 0.577... it the Euler-Mascheroni constant. The amplitude as shown in Fig. [T] 
drops to this level during a time t ~ 1/E g , indicating the time permitted by the energy- 
time uncertainty relation over which virtual excitations above the gap edge may be 
maintained. 
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The same physics holds for soft gaps such as in the standard superohmic case, 
J q (uj) = a q uj q oj]r q e~ u) l U)c . This model describes phonon baths as well as electromagnetic 
environments blocked off at low frequency by a serial capacitance [6]. Here, K(t) can be 
given in closed form [32]. For q > 3 we find a visibility of V[0, 0] = exp [2a q T(q — 1)] for 
k B T <C uj c . For 1 < q < 3, the amplitude decays non-exponentially, also at long times, 
given by 

e -m ^ exp ( 2a q TT(q - 2) sin ( 




and for q = 2 we find e~ K ^ — > e~ 2a2 ^ 1+K - ) (^^y) 2 • These sub-exponential decay 
features cannot be characterized by a time scale T 2 and may resemble reduced visibility 
if the observation does not span enough time. 

Loss of visibility is not restricted to gapped models. It is usually joined by 
genuine exponential T 2 decoherence. This is illustrated in the Ohmic Spin-Boson model 
J\{uS) = aidje~ w / Wc , which describes the quantum version of classical friction and e.g. 
decribes standard resistive environments as well as gapless electron-hole excitations. 
At finite T, we identify a finite 1/T 2 = a{T. The full shape of the envelope in the 
long-time limit, tT 3> 1, is however e~ K ^ — > (T/cu c ) a e~ t ^ T2 , i.e. there is a visibility 
prefactor v[0, 0] = (T/u c ) a . At T = 0, the long-time limit is never reached as there is no 
low-energy scale in the problem, and we find power law decay, e~ K ^ = (1 + w 2 t 2 ) a ^ 2 . 
Thus, at T — > 0, the very small f[0,0] reflects the fact that for most of the time the 
power-law decay dominates and the long-time expansion involving T 2 becomes valid at 
extremely long times only. The finite temperature behavior is illustrated in Fig. |2j Note 
that this initial decay followed by exponential decay has been discussed by an entirely 
different approach in Ref. [33J. 

3.2. Entangled initial condition 

Factorized initial conditions are standard assumptions in the theory of open quantum 
systems, however, they are not always realistic. It is already seen in our above result 



Eq. (12), that our prediction of an initial drop qualitatively holds for essentially any 
nonequilbrium initial condition. As a complementary case, we chose the af such that 



in eq. (12) Ui = —1, Vi = 0. Physically, this state corresponds to a qubit dressed by 
environmental oscillators at all frequencies. It is the variational (in u and v) ground 
state as long as E <^ f dujJ(uj)/uj. It is thus a two-state analog to the initial condition 



discussed in Ref. [34J. Inspection of eq. (12) shows that these initial conditions are 



special in two ways: i) the envelope e~ K ^ is constant in time, 1/T 2 = and ii) 
visibility V[— 1,0] has the maximum possible value. It is related to the visibility of 
the factorized system by V[— 1,0] = \Jv\0, 0]. These observations can be physically 
motivated from the minimization of the total energy, which implies that no further 
rearrangement of the bath is necessary for forming the proper dressed states and the 
dressing is optimum by accomodating the smallest possible total energy. Consequently, 
any other initial condition is not optimal and contains irreversible parts of the interaction 
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Figure 2. Coherent oscillations and their envelope in the Ohmic spin-boson model 
with ui c = 10c, a — 0.1 and fcgT = 0.5e. The envelope predicted by Bloch-Redficld 
theory has the correct slope at long times, but misses overall amplitude. 

with the environment, i.e., genuine T 2 decoherence (i.e. finite T 2 ), which we have already 
identified above as redistribution of surplus energy. 

We appreciate from these results, that the initial conditions play a decisive role. In 
fact, the visibility prefactor is a long-time consequence of short-time physics governed 
by the initial conditions. The choice of initial condition depends on the type of physical 
environment and experimental procedure: A detector is typically switched on during 
the experiment and is well described by factorizing initial conditions. If the switch 
is non-adiabatic and it happens on a time-scale r c , this can be taken into account 
by chosing u(u) = for ujt c <C 1 and u{uo) = 1 for ujt c ^> 1 as an initial state 
in the expression for the visibility. The baths in the material and/or the control 
environment are permanently coupled to the qubit. In a typical experiment, where 
on top of the pure dephasing Hamiltonian eq. Q a large Aa x term is applied at times 
t < 0, the entangled initial condition gives a realistic approximation for the initial state, 
predicting a reduced visibility at all times. Realistic experiments, involving a non- 
trivial preparation sequence encoded in 9(t) at t < may be described by variants of 
these initial conditions. However, the fact that the variational ground state shows the 
highest possible visibility outlines the strength of our result: In many experiments, 
the detector signal representing s x = 1 is not known a priori and is obtained by 
ground state measurements. Following our results, the visibility obtained in a dynamical 
experiment is always smaller, V[u, v] < V[— 1,0]. Thus, our theory does even apply to 
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these data, the loss of visibility is not an artefact of the initial condition, it rather 
is a generic consequence of the fact that quantum computing takes place far from 
thermal equilibrium. In fact, the short time slip accomodates the extra energy of the 
system relative to the dressed ground state. Moreover, experiments with good short- 
time resolution such as the phase qubit setup show an explicit sign of an initial drop of 
the oscillation amplitude [35] . 



4. Beyond pure dephasing 

So far, we have primarily discussed the exactly solvable pure dephasing point of 
Hamiltonian eq. pi), mainly for being able to easily circumvent the Born and Markov 
approximations. Our qualitative results do however not depend on that assumption. 
Although a more general solution may require more elaborate methods such as path 
integral expansions [2T] , flow equations |37] or real-time RG [38J, we can verify 
this conclusion by elementary means for the gapped environment, when uo^ > E. 
As shown before, this environment does not lead to Markovian decoherence and 
thus all decoherence effects are necessarily purely nonmarkovian. In this case, it is 
legitimate to proceed by perturbation theory in E/ui <C 1. The eigenstates of the 
unperturbed, E — 0, Hamiltonian are two-fold degenerate and can be written as 
|V>±,K» = \±)I\ k D {k) (±X l )\n l ) with eigenenergies = E + £(n; + 1/2)^. 

The perturbation Hamiltonian lifts the degeneracy and splits the doubletts into sublevels 

\g, {rit}} = - sin# eff /2|+, {n t }} + cos0 eff /2|- {n,}} (14) 
|e, {rii}) = cos6» eff /2|+, {ri;}) + sin 9 eS /2\-, {rii}). (15) 

Here, the angle depends on the and is given by tan6 ) e fr({^j}) = c({rii}) tan 6* with 
Cj = ({rii}\D(2\i)\{ni}) . This can be viewed as a down-scaling of the effective tunnel 
splitting connected to that specific energies. The shifted energies are E g / e ({rii}) = 
±E eS {{m})/2 + E{{m}) with E eS ({m}) = Ey'cos 2 # + sin 2 ^c({^}) 2 . 

We now use these approximate eigenstates for the Gedanken experiment analogous 
to the pure dephasing case: We prepare the system in the equal superposition \ip ) = 
-^(\g, {0}) + \e, {0})) factorized to the bath and compute the probability of returning 
onto the initial state, S(t) = |(-?/>o !'*/'(£)) | 2 - \i>(t)) is obtained by expanding the l^o) in the 
basis of the approximate eigenstates, \ipo) — J2i(d + ({rii})\e, {rii}) + d-({rii})\g, {rii}). 
The expansion coefficients read 

MRlH^fsin^+cosf^)) (16) 
MM) = (— f^) + - (^)) ' <"> 

where 7 = ({ni}\D(a)\0). This expression nicely illustrates our previous point that the 
initial state is broken up into entangled states that contain significant bath excitations. 
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We now propagate these in time. We find without further approximations 
Sit) = ^ 7 (K})e- J ^ n ^[cos( J E; cff (K})t) + 

{rii} 

sin(£ - 9 eS ) sm(E cS ({ ni })t)] (18) 

We recover the factorizd pure dephasing result in the case of 9 = 0. For the other 
extreme, 9 = ir/2, we recognize that the last term is dropping out but the sum cannot 
be performed analytically. We observe because the {rij} dependence of the argument of 
the cosine provided additional decoherence. For general 9, the second, phase shifted term 
that originates from misalignment of the pure and effective magnetic field 9 ^ 9 e g on 
the Bloch sphere provides yet another dephasing channel. A full numerical analysis 
of these cases goes far beyond the scope of this paper. Qualitatively, we see that 
the dephasing familiar from the pure dephasing case combines with further dephasing 
channels, rendering our previous discussion to remain valid beyond the pure dephasing 
point. 

Note that E e ^ at n« = leads to the same expression that appears in the adiabatic 
renormalization treatment at 9 = 7r/2 [16]. 

Clearly, the results of this section recover the features of the above discussion and 
shows, that the assumption of pure dephasing is not crucial for the physics of the 
loss of visibility. Any more quantitative statement should build on more quantitative 
numerical methods. A wealth of these methods has been developped in the last years, 
including work on the numerical renormalization group [39], analytical RG [3D], and 
non-Markovian master equations [41J. 



5. Possible relevance for superconducting qubits 

We have seen that for experiments with long T2, a pronounced loss of visibility is 
governed by gapped baths, unlike the temporal decay, which is dominated by 1//- 
noise. Such baths can be identified in experiments. We list a few cases and give their 
factorized visibility. Note that higher values of the visibility indicate the inadequacy of 
factorized initial conditions. 

In the flux qubit[7J, the contribution of the junction quasiparticles can be 
approximated by a gapped Ohmic model with cujr = 2A, u c = (R^Cj) -1 where is 
the normal-state resistance of the Josephson junctions and Cj the junction capacitance. 
a = Rq/Rjss ($/$o — 1/2) 2 where -Rq = h/Ae 2 is the quantum resistance for Cooper 
pairs. Using the numbers of Ref. [7] we obtain V[0, 0] = 0.8 

In the Quantronium and other samples, the quasiparticles are largely shunted 
through the capacitor, cu IR > ujq. There, however, an RC-element is fabricated on 
chip |6J. The capacitor is used to decouple the resistor at low frequencies. Indeed, the 
environmental spectral density is super-Ohmic, J{<jj) oc uj 3 at ujRC <C 1. Using the 



expression in Ref. [42J with eq. 12, we obtain a reduction of visibility from this term 
alone of V [0,0] = 0.6. 
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Phase qubits exhibit a large number of spurious resonances. There is no exact 
mapping onto an oscillator bath, however, we can still estimate a visibility logi> = 

IT- 1 2 

— J2i (e-e ) 2 wnere Ei is the energy of the ith resonance and T, is the coupling matrix 
element. This expression can be checked from experimental data when the Tj and Ej 
have been mapped out. 

Recent experiments in charge qubits report extremely high visibility [13]. This 
result agrees with our scenario: The off-resonant degrees of freedom which appear in v 
but not in Ty 2 are filtered out by a cavity. 

The impact of phonons has been studied in Ref. [19]. Due to the super-Ohmic 
spectrum they would be a clear candidate for our scenario, however, the resulting number 
even for factorized initial condition is extremely close to unity. The topic of visibility has 
been studied recently for phase qubits [H] US] and for the Ohmic spin-boson model [22J . 
The connection between exact solutions and approximate master equations is discussed 
in depth in |46j. Both of these work assume weak coupling to the bath and factorizing 
initial conditions. If we expand our results to lowest order, they cover those papers 
as special cases. Similar calculations have been done in Refs. [2U |25j for 9 = and 
factorized initial conditions. 

6. Conclusions 

We have studied the Spin-Boson model beyond the Born and Markov approximations. 
We have shown, that for nonequilibrium initial conditions, a long-time observer will 
measure a reduced visibility of quantum oscillations as a consequence of short-time 
physics involving higher-order processes. This is well compatible with long T 2 for gapped 
or superohmic environmental spectra and with a range of experimental data. 

Discussions with F. Marquardt, S. Kehrein, S. Kohler, D. Vion, J. Clarke, V. 
Golovach, J. von Delft, E. Mucciolo, and L. Hollenberg are gratefully acknowledged. 
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